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Abstract. - We introduce a general framework to study the processivity of molecular motors 
moving along a polar filament and discuss the average time spent attached to the filament 
as a function of a tangentially applied load. Our study of specific models suggests that the 
attachment time of a motor decreases with increasing ATP concentration and that double- 
headed motors such as kinesins lose their processivity under forcing conditions while processive 
single-headed motors are less sensitive to tangential forcing. 



Introduction. - Experiments on motor proteins crucially depend on an important prop- 
erty of such systems, called processivity. Motor proteins such as kinesins, which can hydrolyse 
a large number of ATP molecules before detaching from the microtubule, are called proces- 
sive, whereas those like myosin II, which typically interact just once with actin filaments and 
detach, are called non- processive Q. Physiologically, processivity is important since it per- 
mits transport involving a small number of motors. Non-processive motors, like myosin II 
in muscles, are found to work in large assemblies. After having collected important experi- 
mental information on elementary steps, the velocity as a function of load and the stall force 
of molecular motors such as kinesin ^-Q, there are now efforts to quantitatively study the 
processivity. As a result, there is now ample experimental evidence that applying a force 
tangential to the cytoskeletal filaments and opposing the natural motion of molecular motors 
considerably decreases the processivity . 

A priori, one can define processivity in three different ways. First it can be chosen as 
the average number of chemical cycles before detachment from the filament; second it can be 
defined as the attachment lifetime of the motor to the filament; third, it can be, as in Ref. B, 
chosen to be the mean length spanned by the motor on the filament in a single run. The first 
definition is intrinsic to the process but extremely difficult to measure. The other two are 
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directly accessible to experiment. For the sake of conciseness we focus our discussion on the 
attachment lifetime and only give an example for the run length. 

The existence of a disruptive force normal to the cytoskeletal filament could potentially 
explain the decrease in processivity, since it is clear that tangential forces applied on a bead or 
on a rod attached to the motor produce a normal component on it. The detachment produced 
by normal forces is similar to the detachment of adhesion proteins under load, investigated in 
Rcfs. However, in the case of processive molecular motors there is an additional physical 

mechanism leading to detachment due to applied tangential forces. This effect is the subject 
of the present letter. 

Our calculations show that applying a force strictly tangential to the filament in general 
reduces the average attachment time of molecular motors. This effect is particularly pro- 
nounced for a model motivated by double-headed processive motors, which provides a natural 
explanation for the loss of kinesin processivity observed under loading conditions . It is 
interesting to remark that for a model adapted to single-headed processive motors we find a 
less marked force-dependence of the average attachment time. 

Model. - These results are obtained by taking into account the possibility of detachment 
in iV-state models of molecular motors p-|l6f . In these models, the motor is represented by a 
point particle which can be in one of N internal states (representing the steps of the chemical 
cycle) and is placed at point x in the one-dimensional space representing the filament. For 
each state i, the particle is subject to a periodic potential Wi(x) of period I. The transition 
rates u> ji(x) depend on ATP concentration as discussed in Ref. Jig] . Its hydrolysis drives the 
motion along the filament. 

Detachment can be represented in this picture by allowing the particle to disappear with 
a local detachment rate on{x) which depends both on its state and on its position. It is 
reasonable to expect that a>i(x) ~ exp(Wi(x) / k^T) , since the "attempt rates" for detachment 
should not depend very strongly on position. The exponential dependence is so strong that, 
to a first approximation, one may consider high detachment rates to be concentrated near 
the maxima of the potential; elsewhere we assume a constant, much smaller, detachment 
background. 

We assume that the particle attaches to the filament at t — with a probability distribution 
P^°\x). The behavior of the system can be described by the probability density Pi(x, t) to 
find the particle at position x and state i at time t. This quantity satisfies, for t > 0, the 
evolution equation 



in which fa is the particle mobility and / ext the externally applied tangential force. 

Since the probability that the particle detaches anywhere on the filament between time 
t and t + dt is given by J_°° dx '^2 i ai(x)Pi(x,t), we can evaluate the average attachment 

time r = L°° dtt f_°°dx J2 i ai(x)Pi(x,t). It is straightforward to show that this definition 
is equivalent to 




(1) 



where the currents J; arc defined by 



Ji = -fa [k B Td x P t + (d x Wi - /«») Pi] 



(2) 




(3) 
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Similarly we can define the mean run length as 

/>oo />+oo 

A = / dt dxx^2a l (x)P l (x,t) . (4) 

JO J -oo i 

Introducing the Laplace transform Pi(x,s) = J Q °° dte~ st P(x,t), we obtain 

/+oo 
dxJ2Pi(x,s=0) , (5) 
-oo i 

where Pi{x, s) is, for each s, a solution of the ordinary differential equation 

sPi + ^ + - UijPj) + oaPi = P?\x), (6) 

i 

which we only need to solve for s = 0. 

Method. Since eq. (||) is linear, we can solve it for pj°\x) — 5u S(x — xq) without 
loss of generality. This implies that a motor is initially placed at position xq and in state 
io- In our numerical analysis, we choose xq = and io = 1. Then, for x > xq and x < xq, 
the solution can be found by a transfer matrix technique. The initial solution imposes the 
matching conditions 

Pi(x + e,s) = Pi(xo-e,s) 

Ji(x + e,s) = Ji(x Q - e, s) + 5 iio . (7) 

Assuming that the vector Y = (jjx, y[, . . . , j/jv, Un) contains the values of the solution and its 
first derivatives at point x, one can evaluate the solution at point x + nl as A4 n Y , where the 
transfer matrix A4 is evaluated by explicitly solving eq. (Q) within a period. Notice however 
that, since the boundary conditions imply that P(x, t) vanishes as x — > ±oo at all times, this 
must be also true for P(x,s) for all s. Therefore, when x > xq (x < xq), Y must belong to 
the Af-dimensional space spanned by the eigenvectors of M. {M." 1 ) whose absolute value is 
smaller than 1. The 2N matching conditions (0) are thus in general sufficient to fix the 2N 
coefficients which identify the solution. 

The evaluation of the transfer matrix is difficult in general, but can be made tractable 
by assuming that the tUij(x) and the cti{x) are piecewise constant, and that the Wi{x) are 
piecewise linear. Then the coefficients of eq. ([|) are piecewise constant, and the solution 
can be expressed by exponentials, which must be smoothly connected at the joining points, 
ensuring the continuity of the probability density P and of the current J. Therefore, at the 
price of a special treatment for the period containing xq, the solution of eq. (^|) reduces to 



a set of transcendental equations whose solution can be numerically evaluated jl7,18|. In 
particular, the infinite-force limit |l7]] and the simplest case with only one state can be fully 
solved analytically [la ]. 

Examples. We discuss two simple examples with N = 2 and asymmetric sawtooth 
potentials: (i) Identical shifted states (ISS): This system involves two identical potentials 
shifted by half a period (Wi(x) = Wi{x + 1/2)); (ii) Diffusive steps (DS): In this case we 
choose a flat excited state (W2(x) = const.) and all transition rates for this state are constant. 
This simplified situation is valid as long as the variations of W2 do not exceed kT [0,0. 
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Fig. 1 - Schematic representation of two specific models discussed in the text, (a) Identical shifted 
states (ISS): Transitions are localized in a region of size d^. The detachment rates a are piecewise 
constant with the maximal value a max localized in a region of size d a . Elsewhere a — a m i n . (b) 
Diffusive steps (DS): Excitations are localized at the minima of Wi, deexcitations are homogeneous. 
Detachment rates a are maximal (a ma x) near the potential maximum and a m i n otherwise. The 
potential amplitudes are denoted by U, the asymmetry of the sawtooth potentials is characterized by 
the length a, and I denotes the potential period. The value of W2 can satisfy Wimin < W2 < Wimax, 
but the line has been displaced upwards for graphical clarity. 



The potentials and the transition rates for these situations ISS and DS are defined in FigJfl. 
We denote the potential amplitude by U and characterize the potential asymmetry by a/l, 
where a denotes the position of the maximum of W\ . 

The ISS model is motivated by the motion of double-headed kinesins. Each of the two 
states could correspond to a situation where one of the heads is strongly bound to the filament. 
Because the two kinesin heads are identical, the two states are identical up to a relative shift. 
We choose transitions between the states localized within a region of size d^ near the potential 
minima. At the minimum of W\{x), excitations occur with a rate W21 = Q while 0^12 = lo. 
Because of the symmetry of the states LJ21 = w and ^12 = fi at the minima of W2(x). 
Elsewhere, W12 and LO21 vanish. The rates of detachment from the filament a^x) are chosen 
piecewise constant with a large value ev max located in a region defined by Wi(x) > U — k^,T 
of size d a around the potential maxima, and with a small value a m ; n elsewhere, see Fig.|l|(a). 

Note, that this simple choice captures the essential exponential dependence of the detach- 
ment rate but is more convenient to use for practical calculations. 

The use of the DS model is motivated by the recently observed processive motion of 
artificially constructed individual kinesin heads |L9],|2(]]. The particle, initially lying in the 
potential minima in state 1, gets excited into state 2, where it diffuses. It then falls back 
into state 1, and advances if diffusion has brought it past the nearest potential maximum. 
For this system, we choose excitations that are localized in the potential minimum where 
LJ21 — f2, W21 = elsewhere. The deexcitation rate u>\2 = to is constant along the period. The 
detachment rates oti(x) are large (a\ = a max ) in a region of width d a around the potential 
maxima, and have a small value (ai = tv m in) elsewhere. We choose a.2 to be a constant, 
related to the value of the potential energy in state 2. The fact that the observed attachment 
times of an artificially constructed single-headed kinesin is similar to those characteristic for 
double-headed kinesins [^9| suggests that in this case the potential energy of state 2 is not 
larger than the potential maximum of state 1, and thus that «2 < a max . If the potential 
energy of state 2 is large, one has «2 > a mH , and the corresponding motor is non processive: 
indeed, it is likely to leave the filament as soon as it gets excited. This situation could apply 
therefore, to myosin II. 
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Fig. 2 - Average attachment time to a filament r as a function of the normalized force /ext//s, where 
/ s is a typical value of the stall force. Identical shifted states (ISS) with U /UbT — 20, I = 16 nm, 
a/l = 0.2, /hi = iii = 171 nm pN _1 s~\ w = 10" 4 s" 1 , du/l = 0.02, d a /Z = 0.05, a min = 0.1 s~\ 
amax/amin = 4.8 10 8 for Q = 5.0 10 4 s _1 and £7 = 50 s _1 . For these parameter values, the system has, 
for large Q, a saturated spontaneous average velocity vo — 1.1 /im/s in the absence of external forces 
and the force is normalized with / s = 6.0 pN. Inset: mean run length as a function of the external 
force /cxt, in the ISS with identical parameters. Note the similarity with experimental results of Ref. 
[6]. 

Results. - Figure || shows the average attachment time r as a function of the externally 
applied force for the ISS model for two different values of the excitation rate fl which is an 
increasing function of the ATP concentration. The maximal attachment time is independent 
of SI and occurs in the absence of applied forces. In this regime, the system is processive 
and spends an important amount of time attached. The time r decreases if an external force 
which opposes motion is applied. At stalling conditions, the time of attachment has dropped 
by almost two orders of magnitude: it decreases by several orders of magnitude for larger 
forces. The situation for the DS model is different, as shown in Fig.[|. Forces up to 10-20 
times larger than the stall force are necessary for reducing the bound time by a factor of 10 2 for 
SI = 10 2 s _1 , for SI = 10 5 s _1 the reduction does not exceed a factor of 2 even for large forces. 
The reduced sensitivity of the processivity to external forces in the DS model is due to the 
fact that the x-independent detachment rate «2 governs this regime. Our calculations show 
that the processivity in the absence of external forces decreases for increasing SI (cf. Fig.||(a)). 
This effect occurs in both models but is very significant in the DS model. For the case of 
large SI we find in the DS model a weak dependence of the attachment time as a function of 
tangential load. Figure ||(b) exhibits the behavior of the DS model for the same parameters 
as in Fig.||(a), but with a 2 > ct max (he., W% > W\). The same general trends are visible, but 
the excitation-rate dependence is strongly increased. The processivity and force sensitivity 
are totally lost for high excitation rate SI. 

Discussion and Conclusion. Our study shows that the processivity of a molecular 
motor is influenced by an externally applied force tangential to the filament. This mechanism 
applies to all situations in which a molecule moves along a linear structure. The details of the 
initial condition will not be relevant if the motor undergoes on average several cycles and 
moves a few periods before detaching, like in the case of kinesins. However, for non processive 
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Fig. 3 - (a) Average attachment time to a filament r as a function of the normalized force /ext// s , 
where / s is a typical value of the stall force. Diffusive steps (DS) with U/k-sT = 10, I = 8 nm, 
a/l = 0.2, m = ^ 2 = 5 10 2 nm pN _1 s _1 , u = 100 s _1 , d u /l = 0.02, d a /l = 0.1, a min = 10" 3 s~\ 
amax/amin = 2.2 10 4 , a 2 / a min = 67 (i.e., W 2 ^ 4fe B T), for Q = 1.0 10 5 s" 1 and fi = 1.0 10 2 s -1 . 
In this system, Do — 157 nm/s and / s = 0.33 pN. (b) Average attachment time to a filament r 
as a function of the normalized force / cx t//s, where again / s = 0.33 pN (with the exception of 
the case for £1 — s _1 for which the stall force vanishes). Same model as in Fig.^(a) above, but 
a 2 /a mi n = 1.63 10 5 (i.e., W 2 ~ 12fc B T). 



motors, knowledge of the probability P> is required for making detailed predictions. 

The ISS model suggests that for double-headed motors the time of attachment is particu- 
larly sensitive to the presence of tangential forces. Indeed, in the presence of tangential forces, 
the motor is more likely to be found in high-energy configurations, where the detachment rate 
is large. The DS model demonstrates a different behavior. In this case, the system spends an 
important fraction of time in the second state with a conformation-independent detachment 
rate ai that dominates the detachment probability and sets the order of magnitude of the 
lifetime r. Therefore, the time of attachment depends only weakly on the applied force. 

It has been observed experimentally that the time of attachment of kinesin molecules to 
microtubules decreases as a load is applied tangentially to the filament || . These observation 
are consistent with our calculations corresponding to the ISS model (FigJ^). Equivalently, 
calculations of the mean run length based on the ISS model are fully consistent with the 
results of Ref. ||, (Fig.^j, inset). As a result, kinesins lose their processivity before their 
direction of motion is reversed at stalling conditions and in general reverse motion cannot be 
detected. Furthermore, the attachment time of double-headed kinesins increases by about a 
factor of eight when ATP concentration is decreased at constant load of lpN from 2mMol 
to 5^Mol §. With the choice of Fig. 1(a) we find an increase of the attachment time 
by a factor of 4.5 when decreasing w by a factor that corresponds to the change in ATP 
concentration. If we slightly increase the width of the detachment region a max this factor can 
be further increased to match the observed ratio. Our calculations predict that single-headed 
kinesins, which move processively in the absence of external forces should in contrast to 
double-headed kinesins remain processive in the presence of tangential forces. Therefore, the 
reversal of the direction of motion could be observable under loading conditions. It is further 
gratifying to remark that the DS model is successful in explaining the different processivities 
of KIF1A (Fig^(a)), single-headed conventional kinesin and myosin II (Fig.||(b)): the Wi 
binding energy of KIF1A has to differ by %kgT from that of myosin and by AUbT from that 
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of kinesin. The AksT corresponds indeed to the difference measured in |2Q |. In the myosin II 
case, processivity is totally lost and the calculated behavior in the absence of ATP is fully 
consistent with known experimental data, with and without external force pl|-]23|] . 
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